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Abstract 

Background: Mapping the intracellular fluxes for established mammalian cell lines becomes increasingly important 
for scientific and economic reasons. However, this is being hampered by the high complexity of metabolic 
networks, particularly concerning compartmentation. 

Results: Intracellular fluxes of the CHO-K1 cell line central carbon metabolism were successfully determined for a 
complex network using non-stationary 13 C metabolic flux analysis. Mass isotopomers of extracellular metabolites 
were determined using [U- 13 C 6 ] glucose as labeled substrate. Metabolic compartmentation and extracellular 
transport reversibility proved essential to successfully reproduce the dynamics of the labeling patterns. Alanine 
and pyruvate reversibility changed dynamically even if their net production fluxes remained constant. Cataplerotic 
fluxes of cytosolic phosphoenolpyruvate carboxykinase and mitochondrial malic enzyme and pyruvate carboxylase 
were successfully determined. Glycolytic pyruvate channeling to lactate was modeled by including a separate pyruvate 
pool. In the exponential growth phase, alanine, glycine and glutamate were excreted, and glutamine, aspartate, 
asparagine and serine were taken up; however, all these amino acids except asparagine were exchanged reversibly 
with the media. High fluxes were determined in the pentose phosphate pathway and the TCA cycle. The latter was 
fueled mainly by glucose but also by amino acid catabolism. 

Conclusions: The CHO-K1 central metabolism in controlled batch culture proves to be robust. It has the main 
purpose to ensure fast growth on a mixture of substrates and also to mitigate oxidative stress. It achieves this by 
using compartmentation to control NADPH and NADH availability and by simultaneous synthesis and catabolism 
of amino acids. 

Keywords: Chinese hamster ovary cells, Metabolic flux analysis, Mitochondria, Compartmentation, Mammalian 
metabolism, CHO, Mammalian cell culture, Metabolic transport, Reversibility 



Background carbon sources simultaneously and maintain a stable 

Economic importance and ease of cultivation make metabolism for long periods in batch cultivations. This 
CHO cells a desirable candidate for metabolic studies has led to a wealth of studies aimed at exploring CHO 
in eukaryotic systems. Alongside with being the most metabolism. After the decoding of the CHO-K1 cell line 
important mammalian cell line for producing biopharma- genome [5,6] and transcriptome [7], we can expect such 
ceuticals [1-3], CHO cells are able to grow in suspension studies to increase both in number and complexity, 
cultures using chemically defined media [4], use multiple Metabolic flux analysis has been applied for mammalian 

cells for a long time already but mostly only using metab- 
olite balancing [8]. The labelling of extracellular excreted 
lactate and C0 2 has already been used in the past but only 
at metabolic steady state [8-11]. The use of other labelled 
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metabolites, e.g. amino acids, is more complex because 
of the often high reversible exchange with the media 
[8]. Metabolic flux analysis (MFA) of CHO cell cultures 
evolved from flux balancing analysis [12] to more com- 
plex metabolic or isotopomer dynamic models [13]. 
Newer studies rely on 13 C-MFA applied by fitting the 
summed fractional labeling [14] or by fitting steady-state 
labeling data [15] resulted from using in parallel more 
labeled substrates for determining the intracellular fluxes 
at metabolic steady state in different growth phases. How- 
ever, the labeling patterns of the intracellular metabolites 
or of amino acids from hydrolyzed proteins that are 
usually needed for non-stationary C-MFA are obtained 
through a tedious methodology [16,17] and are susceptible 
to errors stemming mostly from the quenching/extraction 
phase [18,19]. In the absence of metabolite exchange with 
the media, intracellular labeling would reach steady state 
relatively fast, in the order of minutes for glycolytic inter- 
mediates and few hours for TCA cycle metabolites, as it 
was determined in Pichia pastoris [20]. In mammalian 
cells, exchange with the extracellular pools [21] delays the 
intracellular isotopic steady state usually beyond the possi- 
bility to maintain metabolic steady state. Due to the large 
extracellular pools of amino acids, their exchange will 
transfer the time constant of the extracellular labeling 
process, which is in the order of days, to the intracellular 
labeling. One option is to use isotopic non-stationary 
metabolic flux analysis (Inst- CMFA) applied at short 
time scales [22], but this approach has the drawback of 
requiring accurate determination of intracellular concen- 
trations of metabolites [23] . 

Metabolite and reaction compartmentation is important 
for a realistic representation of the mammalian cell 
metabolism, but determining it is more challenging both 
concerning experimental and modeling procedures, as we 
have already reviewed in [24]. In the exponential growth 
phase, a typical culture of CHO is characterized by high 
uptake rates of glucose and glutamine, the Warburg effect 
and the exchange of non-essential amino acids with the 
extracellular media [4,25]. We can expect that by feeding 
a 13 C labeled substrate, some of the extracellular me- 
tabolites will exhibit labeling patterns that can then be 
detected using GC-MS. As these metabolites will be 
enriched in C dynamically, non-stationary C metabolic 
flux analysis (Inst- 13 CMFA) applied to extracellular and 
intracellular isotopomers [26-28] provides a suitable 
framework to determine the intracellular fluxes. Extra- 
cellular pools have a large time scale for labeling (hours) 
compared to the intracellular pools (seconds/minutes), 
thus removing the need to sample intracellular pools 
provided that the labeling information in the extracellular 
metabolites is sufficient. 

We show that by using only the labeling patterns of 
extracellular metabolites produced by feeding [U- 13 C 6 ] 



glucose as the only labeled substrate, intracellular fluxes 
can be successfully determined in a complex, compart- 
mented metabolic network of the CHO-K1 cell line. In 
parallel, we prove that a simplified, non-compartmented 
model is not sufficient for describing the metabolism. 
We also underline the importance of considering reversibil- 
ity when dealing with non-stationary isotopomer models. 

Methods 

Cell culture and experimental set-up 

The CHO-K1 cell line was kindly provided by the Institute 
of Cell Culture Technology (University Bielefeld, AG Noll, 
Germany). The cells were growing in suspension under 
serum and protein free conditions in the chemically de- 
fined medium TC-42 (TeutoCell AG, Bielefeld, Germany) 
supplemented with 4 mM L-glutamine (PAA, Germany). 
Precultures were cultivated in 125 mL baffled Erlenmeyer 
flasks (Corning Inc., Germany) at an initial cell density of 
0.4x10 cells/mL and a working volume of 50 mL on a 
shaking device (225 rpm) at 37°C and 5% CO2 in a humid 
atmosphere. For the main cultivation, cells were harvested 
during the exponential growth phase at a viability of > 98% 
and resuspended in TC-42 medium with 100% [U- 13 Ce] 
glucose (99%, Euriso-Top, Saarbriicken, Germany). The 
main cultivation was performed in a VariolOOO bioreactor 
(Medorex e.K., Norten-Hardenberg, Germany) at batch 
mode with a starting culture volume of 200 mL. The 
bioreactor was inoculated at a cell density of 0.4x10 
cells/ mL. The cultivation temperature was kept constant 
at 37°C and the impeller (3-blade marine propeller) speed 
was set to 300 rpm. During the cultivation, the pH value 
was controlled at 7.2 by gassing with CO2 and by using 
0.5 M sodium carbonate solution. Dissolved oxygen was 
maintained at 30% of the saturation concentration. Sam- 
ples were taken three times a day. Cell density and viabil- 
ity were determined by cell counting using the Trypan 
blue exclusion method. Supernatants were transferred into 
fresh tubes and stored at -20°C until further analysis. The 
average cell diameter was determined using an automated 
cell counter (Invitrogen, Darmstadt, Germany) in a separ- 
ate experiment. This experiment was performed in a shak- 
ing incubator (2 inches orbit, 185 rpm, 37°C, 5% C0 2 
supply) using 250 mL baffled Erlenmeyer flasks (Corning 
Inc., Germany), an initial cell density of 0.4xl0 6 cells/mL, 
a working volume of 100 mL and using the same TC-42 
medium (TeutoCell, Bielefeld, Germany) supplemented 
with 4 mM glutamine. Differences of cell diameters during 
the cultivation were maximum 5% and not taken into 
account. Cell volume was computed assuming the cells are 
spherical using a diameter of 10.6 um. Glutamine deg- 
radation kinetics were determined experimentally in a 
cell-free setup identical to the one employed for cell vol- 
ume estimation. The determined glutamine degradation 
rate constant was kd GLN = 0.0033 IT 1 . 
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Quantification of metabolites 

Quantification of glucose, organic acids and amino acids 
via HPLC was carried out as described previously by 
Strigun et al. [29]. 

Analysis of isotopomer labeling patterns 
Sample preparation 

For determination of labeling patterns of lactate and amino 
acids, 50 ul of supernatants were lyophilized, resolved in 
50 ul N,N-dimethylformamide (0.1% pyridine) and incu- 
bated at 80°C for 30 min. 50 ul N-methyl-N-t-butyldi- 
methylsilyl-trifluoro-acetamide (MBDSTFA) was added 
followed by another incubation at 80°C for 30 min for 
derivatization of metabolites into corresponding dimethyl- 
t-butylsilyl derivatives. For determination of the labeling 
pattern of pyruvate, lyophilized supernatants were resolved 
in 50 ul pyridine containing 20 mg/ml methoxyamine 
hydrochloride and 50 ul MSTFA (Macherey-Nagel, Diiren, 
Deutschland) and incubated at 80°C for 30 min for deriva- 
tization into the methoxyamine -trimethylsilyl derivative. 
Derivatized samples were centrifuged at 13000 x g for 
5 min at 4°C and supernatants transferred into fresh 
glass vials with micro inlets. 

GC-MS measurements 

Extracellular 13 C-labeling dynamics were analyzed by gas 
chromatography mass spectrometry (GC-MS). The GC-MS 
measurements were carried out on a GC (HP 6890, Hewlett 
Packard, Paolo Alto, CA, USA) equipped with an HP5MS 
capillary column (5% phenyl-methyl-siloxane diphenylpoly- 
siloxane, 30 m x 0.25 mm x 0.25 um, Agilent Technologies, 
Waldbronn, Germany), electron impact ionization at 70 eV, 
and a quadrupole detector (Agilent Technologies). The 
injection volume was 1 ul (7683B Autosampler, Agilent, 
Waldbronn, Germany; PTV-Injektor, Gerstel, Miihlheim 
a. d. Ruhr, Germany). Helium was used as carrier gas at 
a flow rate of 1.1 ml/min for analysis of lactate and 
amino acids or 0.7 ml/min for pyruvate analysis. The 
following temperature gradient was applied for lactate 
and amino acid analysis: 135°C for 7 min, 10°C/min up 
to 162°C, 7°C/min up to 170°C, 10°C/min up to 325°C, 
325°C for 2.5 min; inlet temperature: 140°C and heating 
with 720°C/min up to 320°C; interface temperature 320°C; 
quadrupole temperature 150°C The temperature gradient 
for pyruvate analysis was as follows: 70°C for 1 min, 
l°C/min up to 75°C, 5°C/min up to 315°C, 25°C/min up 
to 340°C, 340°C for 5 min; inlet temperature: 70°C 
and heating with 360°C/min up to 360°C; interface 
temperature 320°C; quadrupole temperature 280°C. 

Data analysis 

After identification of metabolites in the scan mode using 
the NIST data bank, quantification of labeling enrichment 
was done in SIM (single ion monitoring) mode in two 



technical replicates using the following unique fragments 
(m/z) containing the complete carbon skeleton of metabo- 
lites: pyruvate 174, lactate 261, alanine 260, glycine 246, 
serine 390, aspartate 418, glutamate 432, glutamine 
431. Mass isotopomer distributions were corrected for 
naturally occurring isotopes using the method of Yang 
et al. [30]. 

Metabolic network models 

Two metabolic networks were established based on ex- 
perimental observations related to metabolite uptake and 
production and extracellular labeling. Both networks 
included: glycolysis; TCA cycle; anaplerotic reactions; 
synthesis of fatty acids, proteins and carbohydrates for 
biomass production; amino acid production and deg- 
radation. Transport from the extracellular media was 
reversible in both models for all metabolites with the 
exception of glucose, asparagine and essential amino acids. 
Mitochondrial transport of malate, oc-ketoglutarate, alanine, 
and reactions of transaminase, malate and lactate de- 
hydrogenase were also reversible. The first model shown 
in Figure 1A considers the intracellular space without 
compartmentation. In the second model (Figure IB), 
the mitochondrial reactions and pools are separated 
from the cytosol. Both models start from the annotation 
of the genomes of CHO-K1 and Mus musculus [31-33]. 
Enzyme localization was established using information 
from the MGI database and data from J. Wahrheit [34] 
who measured compartmented enzyme activity using a 
method adapted from Niklas et al. [35]. Mitochondrial 
transport of alanine was included to explain the existence 
of alanine aminotransferases in both compartments. Me- 
tabolite pools were lumped where it did not influence the 
simulated labeling dynamics. Pentose phosphate pathway 
was reduced to one reaction where one carbon atom is 
lost for each G6P molecule and 5/3 molecules of PG are 
produced. Glycolysis was lumped to three fluxes trans- 
forming G6P into PYR cyt . Isocitrate and citrate were con- 
densed into one citrate (CIT) pool. Succinate, fumarate 
and malate were condensed into one MAL pool. Two 
cytosolic pyruvate pools were used to describe metabolic 
channeling to lactate. Non-essential amino acids catabol- 
ism was lumped to three fluxes fueling the MAL, AcoA 
and GLU pool respectively. No carbon mapping was re- 
quired in this case as essential amino acids are unlabeled. 
Glutaminase activity was mitochondrial [5] and glutamine 
synthetase was cytosolic [36] . Fatty acids, protein and stor- 
age carbohydrates composition of the cell was taken from 
Altamirano et al. [12]. 

In total, the compartmented model consisted of 60 
fluxes and 25 metabolites and the non-compartmented 
model of 42 fluxes and 16 metabolites. The complete 
flux list for the two models, together with the carbon 
transfer rules, is provided in the Additional file 1. 



Nicolae et al. BMC Systems Biology 2014, 8:50 
http://www.biomedcentral.eom/1752-0509/8/50 



Page 4 of 15 




GLN, 



Figure 1 Non-compratmented (A) and compartmented (B) networks of the CHO-K1 central metabolism used for simulations. 

Irreversibility is indicated by simple arrows, and reversibility by double arrows. The reactions depicted in 1B are listed in detail in the Additiona 
file 1 together with fluxes and reversibilities determined. Subscripts meaning: ex - extracellular; c - cytosolic; m - mitochondrial. Abbreviations: 
AA -amino acids; AcoA - acetyl CoA; AKG - alpha-ketoglutarate; ALA - alanine; ASN - asparagine; ASP - aspartate; CIT/ICI - citrate/isocitrate; 
G6P - glucose 6-phosphate; GLC - glucose; GLN - glutamine; GLU - glutamate; GLY - glycine; MAL - malate; OAA - oxaloacetate; PEP - 
phosphoenolpyruvate; PG - phosphoglycerate; PYR - pyruvate; SEP - serine. 



Non-stationary- 13 CMFA methodology 

Isotopic non-stationary metabolic flux analysis (Inst- CMFA) 
comprises: (1) metabolic steady-state balancing of intra- 
cellular metabolites for determining extracellular rates; 
(2) dynamic extracellular metabolite and isotopomer 
balance and (3) dynamic balances of intra-compart- 
mental isotopomers. 

Metabolite balancing 

Net extracellular rates v M ex were determined for each 
extracellular metabolite M_ex for the batch cultivation 
situation, under the assumption of metabolic steady 
state, by fitting the cell density X(t) and extracellular 
concentrations of metabolites C M ex to an exponential 
growth model with specific growth rate ^ (eq. l.a and 
eq. l.b) and constant extracellular rates. Glutamine 
balance included first order degradation in the culture 
media (eq. l.c). 

dX/dt = ftX(t) (l.a) 
dC M _a/dt=v M _ ex X(t) (l.b) 



dC GLN /dt = v GLN -X{t)-kd GLN -C GLN {t) (l.c) 

At intracellular metabolic steady state, the n metabolic 
fluxes that connect the m metabolites are constant and 
satisfy the material balance: 

Gv = 0 (2.a) 

Vj = <Pp j=l..R me as (2.b) 

A, ,=!..„ (2x) 

where G is the m x n stoichiometric matrix and its null 
space v is the vector of net metabolic fluxes which are 
constrained by R meas measured fluxes <pj (2.b) and n 
inequalities (2.c) determined by flux direction. To reduce 
the number of parameters, the free fluxes were extracted 
from the network as described in [37] to produce a de- 
termined stoichiometric system. 

Intracellular and extracellular carbon balance 

The Inst- CMFA framework developed in [23,28] was 
adapted to the case of batch culture cultivation. Isotopo- 
mer balances for extracellular (eq. 3) and intracellular 
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(eq. 4) metabolites were set to be solved together with 
the extracellular balances (eq. l.a, l.b, l.c). 



dlDV M _ ex = 1 

dt Cm ex 



cICm_, 



dt 



•ZDVa 



(3) 



dIDV A 



dt 



J2v r IDV M j n , r vZj n -lDV A 

/=1 



(4) 



with IDV M _ex> IDVjtf cyt, IDV M in , being the isotopomer 
distribution vectors of the extracellular, cytosolic and 
intracompartmental pools of metabolite M. IDV M jn,j is 
the / reaction contribution to isotopomers of metabol- 
ite M, computed using isotopomer mapping matrices as 
described by Schmidt [28]. X conc is the cell volumetric 
concentration expressed in L cell/L media. C M ex is tne 
extracellular concentration of M, is the production 
flux of M expressed in mmol/(L cell x h), v^ c is the 
uptake flux, Vi is one of the R M fluxes entering the 
intracompartmental pool of M, CM_in, and Vffi is the 
flux exiting the pool. The metabolite and isotopomer 
balances from equations 1,2 and 4 are then solved sim- 
ultaneously to obtain the time course of the mass iso- 
topomer distributions. Isotopomer balancing employs 
absolute fluxes that can be computed from the net 
fluxes by introducing a reversibility parameter: 



(5) 



where Vi is the net flux and v,; forward and V; reverse are the 
forward and respectively the reverse exchange fluxes, 

with Vj i forward ^J, reverse ~ ^j, forward — ^ and V '^reverse —0- 

The contribution of reaction / was computed using 
isotopomer mapping matrices [28] that trace carbon 
from the substrate to the reaction products. The initial 
mass distribution of all metabolites was computed consid- 
ering the naturally occurring 13 C fraction (1.1%) and the 
99% atom purity of the employed 13 C labeled substrate. 

The simulated time course of extracellular mass iso- 
topomer distributions (MID) is compared with the experi- 
mental values. The objective function to be minimized was 
expressed as the weighted sum of square differences be- 
tween the experimentally determined and simulated MIDs: 

SSQD = (MIjy im -MID exp ) T -Y^mid ' {MW sim -MID exp ) 

(6) 

where SSQD is the objective function, MID slm is the 
simulated MID, MID exp is the measured MID and E MID 



is the measurement covariance matrix. The optimal solu- 
tion was accepted if it satisfied the x-squared test for 
model verification with 95% probability, and N-p degrees 
of freedom, where N is the number of sampled points (size 
of MID exp ) and p is the number of free parameters. To 
reduce the bias in the objective function generated by very 
small standard deviations, a minimum threshold of 0.005 
was imposed. Accurate confidence intervals and sensitivity 
analysis of fluxes were computed according to [38]. All 
the code was programmed and simulated in Matlab 
[MATLAB and Simulink Release 2013a, The MathWorks, 
Inc., Natick, Massachusetts, United States]. 

Results and discussion 

Metabolic profiling 

The cells exhibited exponential growth for 72 h (Figure 2) 
until glutamine became exhausted and a shift in metabol- 
ism was observed (data not shown). Estimated specific 
growth rate as fitted to eq. (l.a) was fi = 0.0401 h . Up- 
take and production of most metabolites was balanced, i.e. 
they were proportional to growth for the main carbon 
sources and produced metabolites (Figure 2) and for other 
amino acids (Additional file 2). This means that metabolic 
steady state was maintained during the first 72 h of cul- 
tivation. Glucose constituted the main carbon source 
(Table 1), providing 65% of the total carbon entering the 
central carbon metabolism, with an uptake flux of 
371 mmol/(L cell x h). Note that all fluxes are related to 
the cell volume specified by L cell. 39% of the glucose was 
converted to lactate. The observed pyruvate production 
rate was 3.3 mmol/(L cell x h). 

The glutamine uptake flux determined by fitting eq. 
(l.c) to the glutamine concentration over time was 
66.4 mmol/(L cell x h), 16% smaller than that in the 
case when its degradation was ignored. Glutamine uptake 
contributed with 10% to the total carbon pool. The rest of 
the carbon feeding the central carbon metabolism, i.e. 
25%, was obtained from amino acids catabolism. 

Alanine, glycine and glutamate were produced (Figure 3), 
while the other amino acids were taken up in excess of the 
quantity required for biomass synthesis. As a consequence 
of amino acids catabolism, a flux of 92.6 mmol/ (L cell x h) 
fueled the mitochondrial acetyl-CoA pool from the 
degradation of isoleucine, leucine, lysine, phenylalan- 
ine and tyrosine, while a flux of 39.1 mmol/(L cell x h) 
cytosolic AcoA was directed towards fatty acids syn- 
thesis. The catabolism of excess isoleucine, methio- 
nine, phenylalanine, threonine, tyrosine and valine that 
remained after protein synthesis produced 49.6 mmol/ 
(L cell x h) succinate and fumarate (lumped into one 
pool of four carbon di-carboxylic acids, here represented 
as MAL). 

Metabolite dilution by growth was neglected due to 
negligible influence on the total mass balance. 
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Figure 2 Culture profile of the CHO-K1 cells for the first 72 h 
during the exponential growth phase. Experimental values are 
shown with circles and calculated values are represented by 
solid lines. 



Non-stationary labeling experiment 

The MID of extracellular pyruvate, lactate, alanine, 
glutamate, glutamine, aspartate, serine and glycine was 
sampled at 1, 18, 24, 42, 48, 66 and 72 h. During the 
exponential growth phase of 72 h none of the labeling 
reached steady state as shown in Figure 4. Lactate and 
pyruvate exhibited similar labeling dynamics, however 
with different MIDs towards the end of the growth. 
This is surprising since lactate is obtained from pyru- 
vate through the lactate dehydrogenase reaction. The 
predominant lactate fraction, i.e. M + 3, increased to 
0.85 and the pyruvate M + 3 fraction stabilized at 0.81, 



pointing towards glycolytic channeling to lactate achieved 
by the localized cooperation of glycolytic enzymes as 
observed in rapidly proliferating cells [39,40]. From the 
produced amino acids, alanine, also derived from pyruvate, 
had a high M + 3 fraction. Glutamate and glycine M + 2 
fractions increased slowly, with most of the change happen- 
ing in the last 24 h due to the high number of producing 
cells present in the media., Extracellular aspartate, glutam- 
ine and serine were found to be labeled although they ex- 
hibited a net uptake. Glutamine fractional labeling, mostly 
the M + 2 isotopomer, increased sharply at the end of the 
phase, when very little glutamine remained in the media 
and the contribution of secreted glutamine played a large 
role to the labeling state. 

Isotopomer fitting 

Both the non-compartmented and compartmented isoto- 
pomer network models (Figure 1) were fitted to the experi- 
mental mass distributions with the goal of determining 
unknown fluxes and reversibilities. The 7 sampling time 
points of the 8 metabolites produced a total number of 252 
experimental MIDs. 

Convergence to the optimal solution is difficult in 
isotopomer models [41] and the parameter space of the 
objective function is marked by a multitude of local 
minima [42], making gradient-based algorithms unreliable. 
Consequently, we applied a global optimization scheme 
that had the following steps: (1) generate an initial random 
population of (40 x p) parameter sets that satisfy con- 
straints using a simulated annealing-based algorithm, (2) 
submit the population to a SO-generations genetic algo- 
rithm optimization, and (3) refine the best solution using 
a trust region reflective algorithm. Convergence to the 
optimal solution was verified by repeating the optimization 
scheme. One simulation took about 3 s, and the 
optimization procedure required about 40 h on a 2.3 GHz 
QuadCore CPU. All the numerical integration and 
optimization algorithms are found in Matlab toolboxes. 

We had initial difficulties in fitting pyruvate and ala- 
nine labeling dynamics. As it was shown that re- 
versibility greatly affects labeling dynamics [43], we 
assumed that the transport reversibility parameter 
changes in time, even if the net fluxes remain constant. 
The decrease with time of pyruvate transport reversibility 
was mechanistically expressed using a hyperbolic function 

revpYR = ™m™ e > where pyruvate transport reversibility 
reVpYR decreases from a starting value rev pYR . To avoid 
division by zero, a negligible correction factor e was intro- 
duced. Alanine transport into the cell intensifies as extra- 
cellular alanine becomes exponentially more abundant. 
Transport reversibility was expressed in this case as 
rev ALA = a ■ exp(/S • time), where a and fi are parameters to 
be determined. 
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Table 1 Carbon sources for the central metabolism of the CHO-K1 cell line in batch culture during the exponential 
growth phase 



Metabolite 


Target intracellular 
metabolite 


Uptake flux 
[mmol/(L cell x h)] 


Uptake flux 
[Cmmol/(L cell x h)] 


Percentage of the 
total carbon-flux 


Glucose 


Pyruvate 


371.0 


2226.2 


64.9 


Glutamine 


A KG 


66.4 


331.9 


9.7 


Ml* 


AcoA 


92.6 


185.2 


5.4 


AA2* 


Ma late 


49.6 


198.4 


5.8 


AA3* 


A KG 


11.6 


58.2 


1.7 


ASP/ASN 


OAA 


68.1 


272.3 


7.9 


Serine 


Pyruvate, glycine 


48.3 


146.4 


4.3 


TOTAL 






3428.53 


100 



The contribution of AA1, AA2 and AA3 amino acid groups considers only catabolism. Abbreviations: AA1: amino acids catabolized to AcoA (isoleucine, leucine, 
lysine, phenylalanine, tyrosine); AA2: amino acids catabolized to four-carbon dicarboxylic acids (isoleucine, methionine, phenylalanine, threonine, tyrosine, valine); 
AA3: amino acids catabolized to glutamate (arginine, histidine, proline); AcoA: acetyl coenzyme A; AKG: alpha-ketoglutarate; ASN: asparagine; ASP: aspartate; 
OAA: oxaloacetate. 

*Excluding requirements for protein synthesis. 



The 24 free parameters of the non-compartmented 
model (Figure 1A) consisted of 5 fluxes, 18 reversibilities 
and the C0 2 pool. At convergence, the model failed to 
fit the data with the minimized SSQD of 1572, larger 
than x (0.95, 252-24) = 264.2. Pyruvate, lactate, alanine 
and glutamate labeling were fit poorly even when trans- 
port reversibility was variable (Figure 4). In consequence, 
the non-compartmented model was rejected. The low la- 
beling content of pyruvate, alanine and lactate simulated 
with the non-compartmented model is explained by the 
lumping of the cytosolic and mitochondrial pyruvate 
pools. As more than 30% of the carbon feeding the TCA 
cycle is not labeled, it is expected that the cataplerotic 
reactions catalyzed by phosphoenolpyruvate carboxyki- 
nase and malic enzyme will produce a large quantity of 
unlabeled pyruvate, which contradicts the experimental 
observations. 

The compartmented model, consisting of 11 free fluxes 
and other 27 free parameters (reversibilities) depicted in 
Figure IB, fitted the data successfully with the minimized 
SSQD = 249.0 slighdy smaller than x (0.95, 252-38) = 
249.13. The complete experimental and simulated data- 
sets, together with the standard deviation of the measured 
mass isotopomer distributions, are listed in Additional file 
3. The poorer fit of the 66 and 72 h time points for lactate 
and 72 h for alanine can be explained by the metabolic 
shift towards the end of the growth phase. From a param- 
eter fitting point of view, exponential growth will add a 
larger contribution in the objective function to the label- 
ing towards the end of the exponential phase compared to 
the beginning of the experiment because the rates of 13 C 
accumulation in the extracellular media will be much 
larger at high cell densities. This is best evidenced in 
Figure 4 where glutamine, aspartate, serine and glycine 
do not become noticeably labeled until 40 h after the 
introduction of the labeled substrate. 



Metabolic fluxes in the CHO-K1 cell line 

Glucose was converted to PG mostly by bypassing glycoly- 
sis (Figure 5) through the pentose phosphate pathway 
(PPP). The estimated PPP flux was 80% of the total molar 
glucose input flux, a high activity contrasting with results 
obtained by Ahn and Antoniewicz for adherendy growing 
CHO cells [14] but observed for hybridoma [44] and 
cancer cells [45] . A wide range of PPP activities, between 
0-160% of the glucose input flux, was determined for a 
highly-productive CHO line in fed-batch cultivation con- 
ditions at different growth phases [15]. The large quan- 
tities of cytosolic NADPH produced through PPP are used 
to drive fatty acids synthesis and possibly to mitigate oxi- 
dative stress by reducing reactive oxygen species [46-50], 
as it has also been proposed by [15]. Overflow to lactate 
comprised 39% of the pyruvate produced from glycolysis. 
From the rest of the cytosolic pyruvate, 42.7 mmol/(L 
cell x h) were converted to alanine, but most of it was 
transported into the mitochondria and converted to 
AcoA. The channeling flux from PEP to lactate was 
122.7 mmol/(L cell x h), accounting for 42% of the total 
lactate being produced. Low reversibility (Figure 5) meant 
no connection between the two cytosolic pyruvate pools 
PYR cl and PYR c2 existed. However, lactate was produced 
from both cytosolic pyruvate pools, indicating that glyco- 
lytic channeling is not the only lactate source in the cell. 
Our possible explanation is that multi-enzyme complexes 
associated to membrane transporters, as characterized by 
Campanella et el. [51], create a micro-compartmented 
environment in the cytosol. Glycolytic enzymes are partly 
associated and partly soluble, resulting in a mixed response 
in the lactate labeling. 

The carbon flux in the TCA cycle originated mainly 
from AcoA produced from glycolytic pyruvate transported 
into the mitochondria (Figure 5), with significant contribu- 
tions from glutamine and essential amino acids catabolism. 
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Figure 3 Compartmentation of the CHO-K1 metabolism and the fate of extracellular metabolites. Net fluxes are indicated on the gray 
arrows in units of mmol product/ (L cell x h), and reversibility parameter defined as reverse flux/net flux is shown in the square brackets (n.d. = not 
determined). The thickness of the gray arrows is proportional to the forward flux (=reverse flux + net flux), and shown qualitatively for the fluxes with 
variable reversibility. Net flux direction is shown by the black arrow heads. Amino acids catabolism is represented as the sum of the differences 
between amino acid uptake flux and flux required for protein production, reported to the metabolite derived from catabolism. Subscripts meaning: 
ex - extracellular; c - cytosolic; m - mitochondrial. Abbreviations: AA exl - isoleucine, leucine, lysine, phenylalanine, tyrosine catabolized to acetyl-CoA; 
AA eX 2 - isoleucine, methionine, phenylalanine, threonine, tyrosine, valine catabolized to fumarate and succinate; AA ex3 - arginine, histidine, proline 
catabolized to glutamate; AcoA - acetyl CoA; ALA - alanine; ASN - asparagine; ASP - aspartate; CIT/ICI - citrate/isocitrate; GLN - glutamine; 
GLU - glutamate; GLY - glycine; MAL - malate; PYR - pyruvate; SER - serine. 



Such a high activity of the TCA cycle and high connectivity 
with the glycolysis is in contrast with some previous re- 
ports of lower activity during exponential growth phase 
[14,52,53] but similar to [54-56]. The differences can be 
assigned mainly to the use of different cell lines and culti- 
vation conditions like media composition, aeration mode, 
pH control and culture type e.g. suspension or immobi- 
lized, batch or fed-batch. The lower lactate/glucose molar 
ratio of 0.78 reported herein means more pyruvate is 
available for use in the TCA cycle, thus making for a 
more efficient metabolism. Gluconeogenesis was active 
through PEP carboxykinase with 10% of the total flux 
entering the PEP pool, a fact explained qualitatively by the 
presence of M + 2 lactate and pyruvate. In the absence of 



gluconeogenesis, only the M and M + 3 mass isotopomers 
of these metabolites would be present after feeding fully 
labeled glucose. Malic enzyme activity was negligible in 
the cytosol, and this is in agreement with compartmented 
enzyme activity observed by J. Wahrheit [34]. This obser- 
vation reaffirms that PPP was the main source of cytosolic 
NADPH. Mitochondrial malic enzyme was highly active, 
producing one third of the total mitochondrial pyruvate. 
However, a part of the mitochondrial pyruvate was 
recycled back into the TCA cycle via pyruvate carboxylase. 
Mitochondrial malate net transport flux was small and 
reversible. This explained the lack of M + 1 labeling in lac- 
tate, alanine and pyruvate that would have been otherwise 
linked to the M + 2 malate isotopomers that are expected 
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metabolic network specified in Figure 1 A. For the other eight plots (D-K), the compartmented model provided in Figure 1 B was simulated. 



to be obtained in the TCA cycle. As a consequence, the 
M + 2 labeling in these metabolites relies on mitochon- 
drial transport of citrate and on the activity of citrate lyase 
producing cytosolic AcoA and OAA, which is then further 
converted to PEP. 

About one third of the total serine was produced from 
PG, using cytosolic glutamate for transamination. Serine 
was exchanged with the media, thus explaining extracel- 
lular labeling of serine. Half of the serine was not used 
for protein synthesis but was reversibly converted to 
glycine and CI units to sustain the high anabolic 
activity. Glycine was then secreted. The remaining 
excess of serine was converted to pyruvate. Alanine was 



synthesized mainly from cytosolic pyruvate in a highly 
reversible reaction. Connectivity between cytosolic and 
mitochondrial alanine pools and the direction of the 
mitochondrial alanine aminotransferase flux could not 
be determined. However, the transport flux of alanine 
to/from mitochondria was confined between -11.6 to 
13.6 mmol/(L cell x h), i.e. ±25% of the alanine produc- 
tion flux. The flux of 68.1 mmol/(L cell x h) from as- 
paragine and aspartate uptake to oxaloacetate was split 
through aspartate aminotransferases between cytosolic 
and mitochondrial oxaloacetate with a 3/1 ratio, but no 
other details could be inferred due to the low labeling 
level in extracellular aspartate. 
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Isocitrate dehydrogenase (IDH) activity in the cytosol 
could not be reliably determined due to lack of informa- 
tion in directly connected metabolites citrate and AKG, 
but also because it affects the labeling pattern in the 
same way as the mitochondrial isozyme. As a result, the 
flux in the CIT m -CIT c -AKG c -AKG m cycle could not be 
determined. Nevertheless, a net activity of cytosolic GDH 
towards producing the high glutamate flux needed for 



cytosolic transamination reactions implies that AKG is 
either produced in the cytosol by IDH or transported from 
the mitochondria into the cytosol. Mitochondrial glutam- 
ate pool was fed by transporting cytosolic glutamate into 
the mitochondria and by glutamine through GLS activity 
at comparable rates. In the mitochondria, glutamate was 
then converted to AKG and fed into the TCA cycle 
through mitochondrial GDH. In the cytosol, the glutamate 
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produced from AKG in the various transaminase reactions 
was partially converted to glutamine, which was then ex- 
changed with the media, leading to the presence of labeled 
glutamine in the media. In conclusion, simultaneous 
degradation and synthesis pathways for glutamine involve 
glutamine uptake, transport into the mitochondria and 
conversion to glutamate, glutamate dehydrogenation to 
AKG, AKG transport to the cytosol or citrate transport 
and citrate conversion to AKG through cytosolic IDH 
activity, conversion of AKG to cytosolic glutamate, and 
cytosolic glutamine synthesis. 

Transport reversibility 

A very important part in modeling the extracellular 
labeling was considering the reversible exchange between 
the intracellular pools and the extracellular media, a 
phenomenon which affects the dynamics of the label- 
ing process. All sampled extracellular non-essential 
amino acids except asparagine and proline, either pro- 
duced or taken up, were exchanged with the culture 
media (Figure 3). Even if the production flux of alanine 
remained constant throughout the cultivation, the fitting 
remained poor for alanine when considering a constant 
reversibility factor. There, the reversibility was estimated 
to increase with time. The function rev ALA = 0.154 • exp 
(0.0359 • time) was used to compute the forward and re- 
verse exchange fluxes (eq. 5), with both parameters being 
determined with a narrow confidence interval (Additional 
file 1). Time is computed in hours. This successfully 
explained the dynamics of alanine labeling. The time con- 
stant of the reversibility function is a value close to the 
specific growth rate, pointing to the fact that alanine re- 
uptake is correlated to the extracellular concentration. 
Serine secretion flux, as computed with eq. 5, was up to 
35 times higher than the net uptake flux. Glycine re- 
uptake flux was 4.6 times the net production flux. Aspar- 
tate, glutamate and glutamine exchange fluxes were in the 
same order with the net uptake/production flux, as 
expressed by the estimated reversibility parameter values 
of about 1. The confidence intervals for the transport 
reversibility parameters are larger than those for fluxes 
because at high reversibilities the labeling becomes less 
sensitive to small changes in reversibility. 

Pyruvate transport reversibility is described by the func- 
tion revpYR = ftrae^oop where time is specified in hours. The 
hyperbolic function implies that at the beginning of the cul- 
tivation, the intense exchange of pyruvate [57] eliminates 
the difference between the labeling of the intracellular and 
extracellular pools. Pyruvate re-uptake decreases because 
pyruvate concentrations changes slightly (Figure 1) while 
lactate accumulates in the media to reach high concen- 
trations and competes with pyruvate for the monocar- 
boxylate transporters [58,59]. Lactate transport reversibility 



parameter could not be estimated because at the beginning 
of the cultivation there is no lactate present in the 
media that could dilute the intracellular pool and affect 
the labeling dynamics. 

Confidence intervals calculation and sensitivity analysis 

Most of the fluxes depicted in Figure 5 were determined 
with narrow confidence intervals. Interval boundaries 
are not symmetrical due to the non-linear characteristics 
of the mathematical model. Determining both the flux 
and exchange in alternative pathways was not possible in 
the case of high reversibility e.g. for determining the com- 
partmentation of alanine metabolism involving reversible 
transaminase reactions. 

The sensitivity coefficients provided quantitative infor- 
mation about the impact of the measured fluxes on the 
estimated flux values (Figure 6A). Sensitivity analysis 
also evidenced correlations between external fluxes and 
network pathways when [U- Ce] glucose is used. In this 
case, the MIDs of metabolites will depend on the inter- 
play between the multitude of non-labeled carbon 
sources and the glucose feed, as opposed to organisms 
that use only one carbon source [60]. The determination 
of anaplerotic fluxes relied on the supply of four carbon 
metabolites from amino acids catabolism. Changes in 
the glutamine uptake flux (Figure 6A) affected most 
fluxes to a large extent. Nevertheless, the high influence 
was mainly computational, as any increase of the flux 
caused depletion of glutamine at the end of the cultiva- 
tion and dramatically different labeling patterns. Glucose 
uptake flux affected the estimation of the PPP and TCA 
cycle fluxes. Errors in measuring glucose concentration 
over time will propagate in the values of these fluxes, as 
the glucose uptake flux determines the fraction of 13 C 
entering the cell. G6P loses one 13 C through oxPPP, 
therefore estimating the split between glycolysis and 
oxPPP depends highly on determining correctly all carbon 
sources. This explains the high sensitivity of the glycoly- 
sis/oxPPP split to all extracellular fluxes. Unexpected cor- 
relations were observed for the glycine production flux 
that influenced most anaplerotic and aminotransferase 
fluxes. Glycine is produced at the expense of serine, which 
is in turn produced from 3-phosphoglycerate, also con- 
verting glutamate to alpha-ketoglutarate during trans- 
amination, and simultaneously converted to pyruvate, 
thus affecting the availability of both cytosolic pyruvate 
and glutamate. 

Local sensitivity of the SSQD to free parameters com- 
puted as the normalized mean deviation of the objective 
function to variations in the estimated parameters shown 
in Figure 6B evidenced the determinable parameters and 
the redundant parameters. Notoriously difficult to deter- 
mine anaplerotic fluxes PEP carboxykinase and mitochon- 
drial malic enzyme induced a noticeable sensitivity in the 
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Figure 6 Sensitivity analysis to measured fluxes and to parameters. The sensitivity represented in the heat map (A) was computed for the 
compartmented network of CHO-K1 when [U- 13 C6l glucose was used as a labeled substrate, Sens* = (dv ! ,/v*)/(dv m /v'^), where Sens*' is the sensitivity 
of the estimated flux v a re-estimated using values of the measured flux v m at the border of the confidence interval. v m is the average measured flux. 
The bar chart below (B) shows the normalized sensitivity of the objective function (SensSSQD) to the free parameters (fluxes and rev = reversibilities). 
The sensitivity was obtained as a mean value of 100 perturbations of each parameter around the estimated value: SensSSQD = (dSSQD/SSQD*)/(dp i /p*), 
where SSQD* is the optimized value of the objective function (eq. 6), and p* is the estimated value of parameter /'. The rates V; correspond to the rates in the 
network shown in Figure 1A. Abbreviations: subscripts: c: cytosolic, ex: extracellular, m: mitochondrial; ALA: alanine; AcoA: acetyl coenzyme 
A; AKG: alpha-ketoglutarate; ASN: asparagine; ASP: aspartate; CIT: citrate; G6P: glucose 6-phosphate; GLC: glucose; GLN: glutamine; GLU: 
glutamate; GLY: glycine; LAC: lactate; MAL: malate; OAA: oxaloacetate; PEP: phosphoenolpyruvate; PG: phosphoglycerate; PYR: pyruvate; 
SER: serine. 



objective function. The increased network connectivity, 
obtained by coupling ALA or ASP deamination to conver- 
sion of AKG to GLU, contributed to this fact. Oppositely, 
most intracellular reversibilities did not influence the par- 
ameter estimation results. This can be easily inferred from 
the fact that while reaction reversibility affects the dynam- 
ics of intracellular isotopomers, it does not mirror in the 
extracellular labeling apart from the reactions altering the 



carbon backbone. Also, in the situation where high values 
of the reversibilities resulted from estimation, local pertur- 
bations around these values will not influence the MIDs. 

Conclusions 

We have shown that intracellular fluxes of the CHO-K1 
cell line central carbon metabolism in batch culture can be 
determined for a complex network by making use solely of 
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the mass isotopomers of extracellular metabolites resulted 
from feeding [U- 13 Ce] glucose as the only labeled substrate. 
To this end, non-stationary 13 C metabolic flux analysis 
proved an effective tool for unraveling important details of 
the CHO-K1 metabolism. Pathway compartmentation, e.g. 
of anaplerotic reactions and amino acid metabolism had 
to be considered for describing the mass isotopomer dis- 
tribution. We reckon that this fact plays an essential role 
in controlling the availability of NADH and NADPH in 
mitochondria and cytosol, but also in facilitating amino 
acid catabolism. A cancer-like high activity of the pentose 
phosphate pathway produced reducing NADPH partly 
to counteract the oxidative stress generated by the 
mitochondrial respiration and partly to fuel fatty acids 
biosynthesis. Cytosolic pyruvate transport is reversible 
thousand-fold compared to the net production flux, 
indicating that although it is not a carbon source, pyru- 
vate creates an extracellular environment [61] most 
probably by functioning as a balancing system for cyto- 
solic NADH [62]. Considering that metabolite exchange 
with the media played a very important role in determin- 
ing the intracellular fluxes, we expect that future 1 CMFA 
studies of mammalian cells metabolism will include this 
essential aspect. Compartmentation controls the simultan- 
eous degradation and production of non-essential amino 
acids. Most likely, the CHO-K1 cells maintain the expo- 
nential growth phase under batch conditions by using a 
well-connected multi-pool system involving metabolite 
and reaction compartmentation, exchange with the media 
and inter-compartment exchange for controlling the me- 
tabolite and cofactor pools. Further studies on localizing 
enzyme and transporter activity together with sampling 
intra-compartmental concentrations would bring valuable 
contributions at elucidating the function of such cycling 
pathways. Accurate enzyme kinetics and thermodynamics 
[63] in mammalian cells would complement the modeling 
using Inst- 1 CMFA with information about reaction dir- 
ection and reversibility. The knowledge gained through 
Inst- 13 CMFA depicts the CHO-K1 central metabolism 
as a robust, highly interconnected network that ensures 
fast growth and mitigates stress generated by reactive 
oxygen species and the accumulation of lactate in the 
culture media. 

Due to the economic importance of CHO cells, efficient 
production processes leading to high product quality with 
minimum effort are of utmost importance. In-depth 
knowledge about CHO metabolism is expected to provide 
valuable assistance in identifying targets for metabolic 
engineering and guiding the design of feeding strategies 
leading to the development of efficient production pro- 
cesses. Overexpression, silencing or knockout of the spe- 
cific glycolytic enzymes that associate with channeling 
glucose to lactate could either be used to study the control 
of the Warburg effect in cancer cells or for improving 



glucose utilization. Because glutamine is a limiting sub- 
strate, overexpressing glutamine synthetase would enable 
cells to run a more efficient energy metabolism, with 
higher fluxes in the TCA cycle. However, as we have 
shown that compartmentation is important in managing 
metabolites, mitochondrial transporters are likely to 
constitute important targets for genetic modifications. 
Inter-compartmental transport of metabolites is a key 
factor in connecting the cytosol and the mitochondria 
energetically and we reckon that modifying the genetic 
expression of transporters will have significant, perhaps 
surprising effects on the overall metabolism. 

Our proposed methodology of sampling the MID only 
in extracellular metabolites for determining intracellular 
fluxes using Inst- 13 CMFA has the potential of broader 
applications, as it circumvents the need to extract intra- 
cellular metabolites and it is non-invasive to cells. The 
information contained in the extracellular mass isotopo- 
mers has a higher resolution compared to the summed 
fractional labeling used previously in [14]. This is suffi- 
cient for resolving a complex metabolic network when 
more metabolites are produced and/or exchanged with 
the culture media. Therefore, we foresee future applications 
in the study of mammalian metabolism at physiological and 
pathological conditions, especially related to compartmen- 
tation, as reviewed in [64], and oxidative stress, e.g. in can- 
cer, neurodegenerative disorders and ageing. Knowledge 
about the metabolism at the compartment level will be 
essential for identifying therapeutic targets and under- 
standing disease mechanisms. Similarly, the method 
could be applied to other enzymatic systems or pro- 
karyotic cells where an extended metabolite exchange 
with the media occurs. 

Additional files 



Additional file 1: List of reactions in the non-compartmented 
central carbon metabolism of CHO-K1 (Table SI). List of metabolic 
reactions, fluxes and reversibilities in the compartmented central carbon 
metabolism of CHO-K1 (Table S2). Carbon transfer rules are provided in 
the parentheses after each reaction. Reversible reactions are designated 
by double arrows. Reversibility is computed as the ratio between the 
reverse flux and the net flux. 

Additional file 2: Complete culture profile of CHO-K1 during the 
exponential growth phase.The lines represent the fitted concentration 
profiles to the experimental values (dots) and in the boxes are the 
determined extracellular rates [mmol/(L cell x h )] together with the 95% 
confidence intervals. Glutamine uptake was determined by considering a 
spontaneous degradation rate of 0.0033 h _1 .The exponential growth 
phase is shown in the last plot. 

Additional file 3: Experimental and simulated mass isotopomer 
distributions of extracellular metabolites and used standard 
deviations. 
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